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We discuss recent results of the replica approach to statistical mechanics 
of a single classical particle placed in a random A^(^ l)-dimcnsional Gaus- 
sian landscape. The particular attention is paid to the case of landscapes 
with logarithmically growing correlations and to its recent generalisations. 
Those landscapes give rise to a rich multifractal spatial structure of the as- 
sociated Boltzmann-Gibbs measure. We also briefly mention related results 
on counting stationary points of random Gaussian surfaces, as well as on- 
going research on statistical mechanics in a random landscape constructed 
locally by adding many squared Gaussian-distributed terms. 

PAGS numbers: 05.40.-a, 75.10.Nr 

One of the simplest models with quenched disorder - a single classical 
particle subject to a superposition of random Gaussian potential F(x) and 
a non-random confining potential T4on(x), with x G - turns out to be a 
surprisingly rich system, characterised by a non-trivial dynamical behaviour 
as well as interesting thermodynamics. Denoting the total potential energy 
as 'H(x) = ycon(x) -|- V{x), the statistical mechanics of the model is con- 
trolled by the free energy: 



as a function of the inverse temperature /? = 1/T, and the sample size 
L, with brackets standing for the averaging over the Gaussian potential 



* The text is based on presentations at meetings "Random Matrix Theory: from fun- 
damental physics to applications", May 2-6, 2007, Krakow, Poland and "Condensed 
Matter Physics meets High- Energy Physics", July 13-18, 2007, St. Petersburg, Rus- 
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distribution. The covariance function of the random part is usually chosen 
in the form ensuring stationarity and well-defined large- limit: 

(F(xi)F(,X2))^ = 7V/y(^^(xi-X2)=^) . (2) 

Important information about the structure of the Gibbs-Boltzmann equi- 
librium measure ^/^(x) = -^py exp — /?'H(x) can be extracted from the knowl- 
edge of moments 

In the thermodynamic limit of the sample volume V;, ^ oo one expects 
typically 

m, ~ V-^^ , (4) 

where the set of exponents Tq reflects the spatial organization of the Gibbs- 
Boltzmann weights. For example, if the weights are of the same order of 
magnitude across the sample volume, the normalization condition implies 
locally P/3(x) ~ and a simple power counting predicts the exponents 
Tq = q — 1. In such a situation it is conventional to speak about a delo- 
calised measure. The opposite case of a fully localised measure describes the 
situation when essential Gibbs-Boltzmann weights concentrate in the ther- 
modynamic limit in a domain with the finite total volume <C Vl ^ oo, 
and are vanishingly small outside that domain. This situation is obviously 
characterized by trivial exponents Tg>o = and Tg<o = oo. Finally, in 
many interesting situations the exponents Tq may depend on q nonlinearly, 
and in this case one commonly refers to the multifractality of the measure. 
The Eqs.dS]) and ^ imply the following expression for the characteristic 
exponents Tq in the general case 

Tq = \q\PH\Q\l3) - qPHP) (5) 
relating them to the appropriately normalized free energy of the system: 

The investigation of multifractal measures of diverse origin is a very 
active field of research in various branches of physics for more than two 
decades. In recent years important insights were obtained for disorder- 
generated multifractality, see [1] for a comprehensive discussion in the con- 
text of the Anderson localization transition, and also [2] for an example 



krakowO? printed on February 2, 2008 



3 



related to statistical mechanics with disorder. The multifractality of ran- 
dom Gibbs-Boltzmann measures in a context related to ours appeared in 
the insightful paper [3]. A popular way of characterizing multifractality 
invokes the so-called singularity spectrum function f{a). The latter func- 
tion is used to characterize the number dN{a) = V^^^^da of sites in the 
sample where the local Gibbs-Boltzmann measure scales as P/3{r) ~ 
in the thermodynamic limit. The definition allows to extract the typical 
characteristic exponents Tq as, see e.g. [1] 

Performing the a— integration by the Laplace method we obtain that the 
positive values of the multifractality spectrum /(a) are related by the Leg- 
endre transform to the set of exponents r^: 

■Tq = a^q - f{a^), q = f'{a^), (8) 

Thus, the knowledge of the free energy J-{fi) in ([6]) allows one to characterize 
the positive part of the multifractality spectrum of the Boltzmann-Gibbs 
measure. 

Early works by Mezard and Parisi [5] , and Engel [5] used the replica trick 
to calculate the free energy Eq.([T|) of an infinite system, L = oo, confined 
by the simplest parabolic potential Vconi'^) = > 0. Employing the 

so-called Gaussian Variational Ansatz (GVA) the authors revealed the exis- 
tence of a low-temperature phase with broken replica symmetry, hence bro- 
ken ergodicity. They were followed by Franz and Mezard |6j and Cugliandolo 
and Le Doussal [7] papers on the corresponding dynamics revealing long- 
time relaxation, aging, and other effects typical for glassy type of behaviour 
at low enough temperatures. The nature of the low-temperature phase was 
found to be very essentially dependent on the type of correlations in the 
random potential, specified via the covariance function described in Eq.dS]). 
Namely, if the covariance fviu) decayed to zero at large arguments u, the 
description of the low temperature phase was found to require only the 
so-called one-step replica symmetry breaking (IRSB) Parisi pattern. This 
effect correctly captures the statistics of the low-lying minima of associated 
Gaussian energy landscapes [8]. 

In contrast, for the case of long-ranged correlated potentials with fviu) 
growing with x as a power-la\\Q the full infinite-hierarchy Parisi scheme of 



To be more precise, at large separations we require the structure function 
([y (xi) — V {x2)]^) y oc (xi — X2)^^ , < 7 < 1. However one can easily satisfy 
oneself that in the present model under consideration the difference between the 
structure function and the covariance is immaterial for the free-energy calculations. 
This will be no longer the case for the model discussed in the end of this article. 
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replica symmetry breaking (FRSB) had to be used instead. 

Based on formal analogies with the Hartree-Fock method Mezard and 
Parisi[l] argued that GVA-based calculations should become exact in the 
limit of infinite spatial dimension A^. In a recent paper [9j the replicated 
problem was reconsidered in much detail by an alternative method which 
directly exposed the degrees of freedom relevant in the limit N ^ oo, and 
in this way allowed to employ the Laplace (a.k.a. saddle-point) evaluation 
of the integrals. The results obtained in [9j by this method for the parabolic 
confinement case fully reproduced those obtained by GVA in [H [5] . 

The method of [9] also works for statistical mechanics of a single par- 
ticle inside any spherical sample |x| < L of a given radius L which makes 
it particularly suitable for studying, e.g., multifractality of the associated 
Gibbs-Boltzmann measure. To this end it is easy to understand that the ra- 
dius L must be scaled with the dimension as L = R^/N to ensure nontrivial 
results when N ^ oo. The effective size R < oo (which is actually half of 
the length of an edge of the cube inscribed in this sphere) can be used as an 
additional control parameter of the model. In particular, the chosen scaling 
L = RsfN ensures that the sample volume Vl = vr^/^ r(A72+i) retains in 
the limit N ^ oo the natural scaling with size R and dimension A^. Namely, 
for R ^ 1 we have InV^ = A^lni?+ smaller terms, which is very essential 
for the analysis of multifractality. 

In what follows we thus concentrate on the case of no confinement po- 
tential Vconi^) = 0, and choose R to have any fixed value (eventually we 
will be interested in a kind of thermodynamic limit R oo). One of the 
observations made in [9] is the existence of a simple mathematical criterion 
which formally differentiates between the short-range correlated potentials 
and their long-ranged counterparts. Namely, assume the covariance func- 
tion fv{u) in ([2]) to satisfy technical conditions fy{u) < 0, fy{u) > and 
fy{u) < for all u > 0, and also fy{u) — > when u ^ oo. The criterion is 
based on considering a combination A{u) expressed in terms of f{u) a^ 

where dashes indicate the order of derivatives taken. Then any potential 
satisfying A(u) > 0, V« > (this family includes, e.g., the potentials with 
fv{u) = exp [—(a + bu)°'], such that a > 0,6 > and < a < 1) turned 
out to have the low-temperature phase which is necessarily of IRSB type. 
The standard replica stability analysis of this IRSB low-temperature phase 
revealed that the stability is controlled by two eigenmodes, denoted in [U] 



^ This eventually coincides, up to the overall sign, with the standard definition of the 
so-called Schwarzian derivative {f'{u),u}. 
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as Aq and ( see equations (B.29) and (B.30) of the Appendix B of that 
paper). If both are positive, ah other eigenvalues of the stabihty matrix 
are positive and the IRSB solution corresponds to an extremum of the free 
energy functional stable with respect to small variations. And those two 
eigenvalues were indeed found to be strictly positive as long as A{u) > 0. 

The situation was found to be very different for the potentials with 
A{u) < 0, Vii > (this family includes, most notably, the powerlaw-correlated 
potentials with the covariance of the form fv{u) = f i^) — d'^ {u + a)^ , /(O) > 
g'^ , < ^ < 1). The low temperature phase is now of FRSB type, and it 
is only "marginally stable". Indeed, the stability matrix for this type of the 
replica symmetry breaking can be shown to contain always a family of zero 
"replicon" modes, see e.g. jlOj for a calculation in the framework of GVA. 

Finally, the above criterion naturally singles out the random potentials 
satisfying A{u) = 0, Vu > as a boundary case between the two regimes. 
Denoting f{x) = fy{u) and noticing that A{u) = implies f'/f^^"^ = const, 
we find the function fv{u) to be equal to fv{u) = Cq'^ In {Cqu + Ci) + 
C2U + C3, where Cj are arbitrary constants. The condition fy{u) — > when 
u — > 00 then selects the case of logarithmic correlations as the only possible, 
which we write as 

fv{u) = h-gHniu + a"). (10) 

This latter choice turns out to be in many respects the most interesting 
situation. Indeed, in [9] it was shown that it leads to a phase diagram 
which combines features typical for the short-ranged behaviour with others 
characteristic of the long-ranged disorder. As a particular interesting feature 
we would like to mention that although the low-temperature phase can be 
thought of as described by a special case of IRSB breaking scheme, the 
relevant eigenvalues Aq and h*^ of the stability matrix identically vanish 
everywhere in the low-temperature phas^, rendering IRSB phase in this 
special case marginally stable. 

The qualitative difference between the three cases - short-ranged, long- 
ranged, and logarithmic, is most clearly seen in the thermodynamic limit 
of large sample size i? — > 00. One finds that for a typical short-ranged po- 
tential the domain of existence of IRSB phase vanishes as long as ii —> 00. 
For example, the transition (de-Almeida-Thouless[llJ, AT) temperature sig- 
nalling of instability of the replica-symmetric solution typically behaves as 
Tat{R) ~ R^y/f^AW) and rapidly tends to zero for decaying correlations. 
For any fixed temperature T > in the limit i? — > 00 the system is effec- 
tively in the high-temperature replica symmetric phase, and the free energy 

^ This fact, though not exphcitly mentioned in [9], immediately follows from definitions 
(B.29) and (B.30) after substituting for qi — qo = Q and qd — qi = y the expressions 
(74) and (79) of that paper. 
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behaves asymptotically like F{T) ~ —TNlnR. Actually this result can 
be seen as a purely entropic contribution, and in particularly implies via 
Eqs. (|5|6p the trivial scaling of the exponents Tg = q — 1, corresponding to 
the totally delocalized Boltzmann-Gibbs measure. 

In contrast, for a power-law growth of correlations one finds that the low- 
temperature glassy phase occupies bigger and bigger portion of the phase 
diagram with growing radius R. Indeed, the transition temperature can 
be shown to grow with R as Tat{R) ~ R'^ , and increasing the system 
size ill at any fixed temperature T > results in the free energy given 
asymptotically by the temperature-independent value F(T)|/j^oo ~ —NR'^. 
This expression actually coincides with the typical minimum of the energy 
function for our system. The corresponding exponents rg>o = 0. In a sense 
the system of this type is always "frozen" in the thermodynamic limit, and 
indeed the Boltzmann-Gibbs measure is localized on a few deep minima. 

Only for the logarithmic case Eq.Q the transition temperature tends in 
the thermodynamic limit to a finite value Tat{R oo) = g , and the free 
energy asymptotics depends non-trivially on the temperature: 



This is natural to interpret as a freezing transition, precisely of the same sort 
as appeared in the celebrated Random Energy Model (REM) by Derrida[14]. 
The same expression for the free energy appeared actually in studies of a 
zero-energy wavefunction for Dirac particles in dimension N = 2 and ran- 
dom magnetic field [3|, after a mapping to a problem of statistical mechanics. 
The Boltzmann-Gibbs measure in this particular case is characterized via 
a set of non-trivial multifractality exponents Tg, see [3] and also |12j . The 
ensuing multifractality spectrum f{a) is simple parabolic for all tempera- 
tures, and shows interesting "freezing" behaviour for a ^ at T = g, i.e. 
at the point of ergodicity breaking. 

We thus see that our results have counterparts in the finite-dimensional 
systems. Actually, understanding the generic statistical-mechanical be- 
haviour of disordered systems for finite remains very challenging problem. 
To this end, rather detailed attempt of investigating our model for finite di- 
mensions A < oo in the thermodynamic limit L — > oo was undertaken in a 
very insightful paper by Carpentier and Le Doussal p!5]. That paper also 
can be warmly recommended for describing the present model in a broad 
physical context and elucidating its relevance for quite a few other interest- 
ing and important physical systems, as e.g. directed polymers on trees [13] . 
The work was based on employing a kind of real-space renormalisation group 
(RG) treatment augmented with numerical simulations. The authors con- 
cluded that for finite spatial dimensions neither models with short-ranged. 




T>g 
T<g 
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nor with long-ranged correlations can display a true phase transition at 
finite temperatures T > 0. And only if correlations grow logarithmically 
with the distance, for such marginal situation the true REM-like freezing 
transition indeed happens at some finite T > at any dimension > 1. 
Fortunately, the logarithmic growth is not at all an academic oddity, but 
rather emerges in quite a few systems of actual physical interest, see jl2] 
for a detailed discussion and further references. 

We thus see that the picture following from the results of [9j for the 
thermodynamic limit (understood as R ^ oo) of the model in infinite di- 
mension is in overall qualitative agreement with the N < oo renormalization 
group studies of the same model in the limit L ^ oo. Another fact which is 
perhaps worth mentioning is that a recent work [15] claimed that IRSB low- 
temperature phase fails to survive in finite spatial dimensions, the fact being 
related to absence of marginally stable modes in the fluctuation spectrum. 
If one assumes that the validity of that claim extends beyond the particular 
model considered in [15], then in our case IRSB phase in finite dimensions 
has no chance of survival for any short-range potentials, but in the loga- 
rithmic case it could survive due to the mentioned marginal stability. This 
picture would be indeed in agreement with the above-discussed RG results 
of [12]. We consider further work in this direction highly desirable, although 
it is clear that performing any perturbative expansion around N = oo limit 
is expected to be a rather technically challenging task. 

We end up our presentation by giving a brief overview of a few most re- 
cent advances in understanding the statistical mechanics of a single particle 
in random high-dimensional potentials. 



As revealed by J. -P. Bouchaud and the present author in |16^ I17j . the 
picture of potentials with short-ranged, long-ranged, and logarithmic corre- 
lations presented above is still incomplete, and misses a rich class of possible 
behavior that survives in the thermodynamic limit R ^ oo. Namely, given 
any increasing positive function $(y) for < y < 1, one can consider poten- 
tial correlation functions /y(n) in the right-hand side of Eq.([2]) which take 
the following scaling form 



This type of potential can be constructed by a superposition of several 
logarithmically correlated potentials of the type (jlOp with different cutoff 
scales tti, and allowing those cutoff scales to depend on the system size R 
in a power-law way: ~ R'^% < Ui < 1 |16j . 



0.1. Multiscale logarithmic potential 




0<u< R^, 



(12) 
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The thermodynamics of such system in the hmit i? — > oo turns out to 
be precisely equivalent [16] to that of the celebrated Derrida's Generalized 
Random Energy Model (GREM)[18l [19]. The REM-like case Eq.([l0D turns 
out to be only a (rather marginal) representative of this class: ^{y) = g^y. 

The leading term in the equilibrium free energy turns out to be of the 
form F{T) = N In RJ^^T), where for < T < Tat = a/^'(1) 

- HT) = T^T) + '^^^^^ + 2 y^)dy , (13) 

where the parameter u^, is related to the temperature T via the equation 

t2 = $'(i/,) . (14) 
For T > Ta_t the free energy is instead given by 

-r(T)=T+m_m. (15) 

These expressions for the free energy can be given a clear interpretation 
as describing a continuous sequence of "freezing transitions" of REM type, 
with freezing happening on smaller and smaller spatial scales [TOj 117], 

The form Eq. (jl3ll5p can give rise to a rather rich multifractal behaviour 
of the Boltzmann-Gibbs measure [T7]. The associated singularity spectrum 
/(a) calculated via Eq.® is positive in an interval a G {amin, otmax)-, where 
the zeroes amin-,Cimax of the function /(a) are given by 

a™„ = -mP) - W C V^dy , amax = -dHP) + 2(3 [' dy . 

Jo Jo 

The singularity spectrum is symmetric with respect to the midpoint of the 
interval of interest, Um = {oimin + cx.max)/2 = —I3J^{(3) > 0, where it has 
the maximum f{am) = 1 as expected. Close to this maximum, namely, in 
the subinterval a G (a_,a+) with a± = Om i 2A{I3)t^^^ where A{j3) = 
/3^(<I>(1) — $(0)) the singularity spectrum has a simple parabolic shape: 

f[a) = 1- ^^^{a- Umf, a-<a<a+. (16) 

In particular, at the boundaries /(a±) = 1 — (3\rp (<I>(1) — <I'(0)). Note that 
in the REM-like limit <&(y) = g^y we have ajnin/max ~^ '^~/+ the 
parabolic behaviour is the only surviving, in agreement with the results of 

mm- 

At the same time outside the interval of parabolicity the general GREM- 
like model shows a much richer multifractal structure manifesting itself via 
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a quite unusual behaviour of the singularity spectrum close to the zeros 
ctmin, C(max- To illustrate this fact, we consider a broad class of functions 
<I>(y) behaving at small arguments y ^ 1 as ^{y) ~ C'^ y'^^^^ with s > 
and the coefficient < C < oo. In particular, in the limiting case s ^ we 
are back to the old REM-like model. Now we can extract the behaviour of 
the /(a) when approaching the endpoints amin or Omax- It is given by 



/(«) ^ -^J^ - a™.v--l^ , (17) 



where 



Or 



^""^ {pcV2FTi)~^ . (18) 



(s + l)(2s + l 



We see that for any s > the derivative of the singularity spectrum diverges 

— ^ 

as / (a) ~ \o!—amin/max\ "'^'^ ~^ This is very different from the Standard 
behaviour observed in other disordered systems [Tj'2]: /'(a) < oo at zeros of 
f{a). At the level of multifractal exponents Tq this feature is translated to a 

rather unusual behaviour for large enough \q\, namely: Tq — qamin = —acq~~^ 
for q > T/Tat, and a similar formula for q < —T/Tat- Note, that in the 
standard situation one always observes linear behaviour Tq = qamin,max 
starting from some value of \q\, see the formula (2.42) in [Ij and discussions 
around it. 



0.2. Extrema of random landscapes and ergodicity breaking 

Another set of recent works on the random Gaussian model with cor- 
relations specified by Eq.Q which deserves mentioning is a continuing 
attempt |20j to relate the phenomenon of ergodicity breaking occuring at 
the level of statistical mechanics to statistical properties of the minima 
(and other stationary points) of high-dimensional Gaussian random surfaces 
'H(x), see |21) for introduction to the problematic. The authors managed 
to show that for a generic smooth, concave confining potentials V^on(x) 
the condition of the zero-temperature replica symmetry breaking coincides 
with one signalling that both mean total number of stationary points in 
the energy landscape, and the mean number of minima are exponential in 
A^. For a generic system of this sort the (annealed) complexity of minima 
vanishes cubically when approaching the transition, whereas the cumula- 
tive annealed complexity vanishes quadratically. One also can investigate 
how the complexity depends on the index of stationary points [22| 120] . In 
particular, in the vicinity of the transition the saddle-points with a positive 
annealed complexity must be close to minima, as they were found to have 
a vanishing fraction of negative eigenvalues in the corresponding Hessian. 
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0.3. Statistical mechanics in a sum of squared Gaussian- distributed 

potentials 

Finally, let us mention recent work [23j on the statistical mechanics in 
the energy landscape given by 'H(x) = ^x^ + X^j=i W^f(x). Here Wi(x), 
with i = 1, . . . ,K are assumed to be independent, identically distributed 

Gaussian tacfons with .ere mean, the vaoance ^W}) = . and the 



.2 

structure function {\Wi (xi) — Wi (x2)]^|) = 24>w { ^^^~^^'' ^- Note impor- 
tant differences from the Gaussian case: (i) the absence of the factor N in 
front of the (co) variance, in contrast to Eq.([2|), and (ii) necessity of specify- 
ing both functions a and (pw-, as the phase diagram will actually depend on 
both of them, in contrast to the the discussion in the footnote 1. The free 
energy of such a system turns out to have a well-defined large— iV limit pro- 
vided we scale K = Nk, and consider the parameter < k < cxd. Naively 
one may think that the central limit theorem (CLT) would imply that the 
sum of K = 0{N) random terms effectively behaves as a Gaussian po- 
tential. A thorough consideration shows that such a reasoning is however 
deficient for the statistical mechanics problem in hand. Indeed, with low- 
ering the temperature deep minima of the resulting potential start playing 
most prominent role, and the description of those minima goes beyond the 
applicability of CLT. This fact suggests that the statistical mechanics of 
such model may have features rather different from the former Gaussian 
case due to different statistics of deep minima[8j. The dynamics in this case 
may also be rather different, see interesting related results in |24| . 

The free energy can be evaluated in the limit — > oo by extending the 
methods of [9], and the system shows both similarities and dissimilarities 
to the Gaussian case. In particular, the difference between the short-range 
and long-range potentials remains to be important, but manifests itself in 
a somewhat different way. One again finds that the replica-symmetric solu- 
tion is unstable at low enough temperatures. Let us discuss here only the 
simplest case of a short-ranged potential with position-independent variance 

a (^jf^ = a = /vk(0), where fwiu) stands for the covariance function of the 

field W, related to the structure function as (j)w{u) = fwi^) — fw{u). The 
equation for the transition (de-Almeida-Thouless) line T^t(m) is then given 
in terms of the structure function by: 



Y _j_ 't>w{TAT) 



2 - , TAT = TAT/^^ (19) 

it; 

Ta-i 



By investigating this expression one finds that the phase with broken replica 
symmetry may exist only as long as the parameter k exceeds some critical 
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value 



1 



(20) 



K > K, 



1 + 



Moreover, for every such k the curvature of the confining potential must 
satisfy the inequality 



In the case of long-range potentials one has to take into account the 
fact of position-dependent variance, which makes the analysis more compli- 
cated, and the corresponding phase diagram quite intricate. These features 
are currently under investigation [23j. 
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